Electromagnetic energy within magnetic spheres 
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Consider that an incident plane wave is scattered by a homogeneous and isotropic magnetic sphere 
of finite radius. We determine, by means of the rigorous Mie theory, an exact expression for the time- 
averaged electromagnetic energy within this particle. For magnetic scatterers, we find that the value 
of the average internal energy in the resonance picks is much larger than the one associated with 
a scatterer with the same nonmagnetic medium properties. This result is valid even, and specially, 
for low size parameter values. Expressions for the contributions of the radial and angular field 
components to the internal energy are determined. For the analytical study of the weak absorption 
regime, we derive an exact expression for the absorption cross section in terms of the magnetic 
Mie internal coefficients. We stress that although the electromagnetic scattering by particles is a 
well-documented topic, almost no attention has been devoted to magnetic scatterers. Our aim is to 
provide some new analytical results, which can be used for magnetic particles, and emphasize the 
unusual properties of the magnetic scatters, which could be important in some applications. 



I. INTRODUCTION 

The research in magneto-optics, both theoretical and 
experimental, has been mainly devoted to the study of 
magnetic properties of thin films. Magneto-optical effects 
are characterized by the change in the state of light po- 
larization in the presence of magnetic materials, both in 
transmission (Faraday effect) and reflection (Kerr effect). 
Brillouin light scattering technique allows the investiga- 
tion of spin waves in magnetic films and layered struc- 
tures through the light scattering by magnons. Here, 
we are concerned with another feature in the magneto- 
optics research: the electromagnetic (EM) scattering by 
magnetic particles [H-Q- Although the EM scattering 
by particles is a well-documented topic [6MToi| . little at- 
tention has been given to the case of EM scattering by 
magnetic particles. Recently, it has been a growing in- 
terest on photonic band gap materials (PBGs) made of 
ferromagnetic materials, like soft ferrites, at microwave or 
radio frequencies [Ill4l4j . Other important applications 
involving magnetic materials, such as microwave filters, 
metamaterials, high density magnetic recording media, 
have been reported [H, EH ■ nere > the approach we follow 
is the classical one for single Mie scattering [f| 0, H, EH > 
in which no applied external field is considered. 

The EM radiation scattering by magnetic spheres is 
described on the basis of the Mie theory, in which an in- 
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cident plane wave, with wavenumber k, is scattered by a 
homogeneous sphere of radius a. We assume that both 
the scatterer and the medium are nonmagneto-optical ac- 
tive and that the incident radiation is a vetorial wave. In 
general, the bulk of analysis takes place in far-field ap- 
proximation, ignoring the evanescence and the internal 
fields in the scattering center The interest, there- 

fore, lies in the behavior of the scattered fields and all the 
quantities of interest to describe EM scattering by spher- 
ical particles, such as cross sections and the anisotropy 
factor (cos 9) (i.e., the mean value of the cosine of the 
scattering angle 9) , can be expressed in terms of the Mie 
coefficients a n and b n Q. For magnetic scatterers, in 
particular, a n and b n have been obtained by Kerker et 
al. [1| . Nevertheless, here, as in the original work of Bott 
and Zdunkowski 17[ for nonmagnetic spheres, the inter- 
nal fields in Mie single scattering gain special attention 
and some related quantities are studied. 

Bott and Zdunkowski [ItJ present the exact and ap- 
proximate analytical expressions for the time-averaged 
EM energy within a dielectric sphere. The calculations 
have been anchored on the rigorous Mie theory, and the 
expressions have been derived, as usual, with the assump- 
tion of equality between the magnetic permeability ten- 
sors of medium and particle. This configuration is de- 
nominated nonmagnetic scattering It is pointed out 
111 that those calculations are of importance for the 
study of photochemical reactions within atmospheric wa- 
ter spheres. 

The aim of this paper is to provide a detailed descrip- 
tion of the time-averaged EM energy within magnetic 
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particles (assumed to be spherical), emphasizing their 
unusual properties, which in turn could be explored in as 
microwave filters and PBGs [15] or in the search ofpho- 
ton localization in the multiple scattering regime 0, @ ■ 
In Sec. |TT] is presented a brief resume about the construc- 
tion of the exact solution in the single magnetic Mie scat- 
tering and its principal analytical results. Both the EM 
internal fields and the magnetic Mie coefficients are pre- 
sented. The determination of the exact expression for the 
time-averaged EM energy within a magnetic scatterer is 
shown in Sec. IIII1 The problem symmetry allows us to 
express separately the contribution of the radial and an- 
gular components to the average internal energy. A new 
expression for the absorption cross section in the mag- 
netic case is determined. To validate our expressions, 
for instance, we determine the same particular relations 
studied in [17| . Special attention is paid to our approach 
concerning the differences to [13 . Finally, we present 
some numerical results in Sec. lIVI We compare the mag- 
netic and nonmagnetic scattering. The basic relations 
involving the Bessel and associated Legendre functions 
are presented in the Appendix [X] Those expressions are 
important to the calculation of the quantities related to 
the time-averaged internal energy. In the Appendix [B] 
some classical limiting cases are considered and we give 
a set of approximated magnetic Mie coefficients. 



II. ANALYTICAL CALCULATION OF 
SCATTERING QUANTITIES 

To deal with EM wave scattering by a single particle 
embedded in a medium, one must assume some special 
features for the medium and the incident wave. Among 
these assumptions, the particle is considered isolated in 
an infinite medium, which allows one to ignore the effect 
of multiple scattering 0, 0] ■ Both particle and medium 
are considered linear, homogeneous, and isotropic, hav- 
ing inductive capacities (ei,/ii) and (e,/i), respectively. 
Thereby, once we assume the media are nonmagneto- 
optic active, those tensors, respective to magnetic per- 
meability (/i) and electric permittivity (e), can be ex- 
pressed by a scalar quantity times an unitary tensor. In 
particular, it is assumed that there are absorptive com- 
ponents within the scatterer, so the quantities ei and /ii 
are complex. 

The incident radiation is considered plane, monochro- 
matic, and polarized complex EM wave, which is ex- 
pressed as 

Ei(r,i) = E exp[z(k-r-orf)], (1) 

with wave amplitude Eo = Eoe x , wave vector k = ke z 
and angular frequency to = 2ir/\ 7 where A is the wave- 
length. Due to the spherical symmetry of the scattering 
center, there is no loss of generality taking the electric 
field polarized on the x axis direction. Also, the linearity 
of the macroscopic Maxwell equations and Fourier the- 



ory allow one to generalize this monochromatic case to a 
polychromatic one @. 

The incident, scattered and internal vector waves have 
the same angular frequency w, once we are not account- 
ing for possible energy variations in the interaction with 
the scatterer. Thus, quantum fluctuations such as in Ra- 
man scattering is neglected, and a classical description is 
adopted 

In the rigorous Mie theory it is quite common to as- 
sume the equality between the magnetic permeability 
tensors of the particle and medium. This consideration 
ignores the most general case in which these complex 
tensors are different. The absolute value of the magnetic 
permeability \x\ can assume values much larger than /z, 
as in the case of soft ferroma gne tic particles in the mi- 
crowave range, for instance [2I In this present work, 
the Mie coefficients are recalculated in this general case, 
referred to as magnetic scattering 0, , and the asso- 
ciated Mie coefficients of the internal fields, which have 
not been studied so far, are presented. The expressions 
here obtained are valid for a wide class of soft ferrites 
with magnetic loss. The assumption of the isotropic mag- 
netic permeability (and electric permittivity) allows one 
to solve the scattering problem in a simple way. However, 
to lower the magnetic loss of these magnetic materials, 
it is usually to consider them in presence of an applied 
external magnetic field [IlTfl3l [l5T[l6j. In this situation, 
the relative magnetic permeability is anisotropic and its 
tensor elements are depends on this externally applied 
field. 

From the Maxwell theory, we have that a time- 
harmonic EM field (E, H) in an homogeneous, isotropic 
and linear medium must satisfy the vectorial Helmholtz 
equation [V 2 + fc 2 ]E = 0, [V 2 + fc 2 ]H = 0, where 
k 2 = — e/j,d 2 = w 2 e/it, and be divergence-free null: V • E = 
V • H = 0. The quantity |k| = k is the wavenumber and 
it is related to the travelling wave. In addition, E and H 
are not independent: VxE = zw/iH, V x H = — jweE. 
To simplify the resolution of these equations, we build so- 
lutions which are dependent of a scalar function called 
a generating function for the vector harmonics 
this particular case, once the symmetry of the problem 
is spherical, the solutions of the equations above are the 
spherical vector harmonics, expressed by M = V x (rip) 
and N = V x M/fc. The imposition of the vector har- 
monics are solutions of the Maxwell equations implies 
that [V 2 + k 2 ]ip = 0. Thus, the problem of the scattered 
waves by a spherical particle resumes to solve this scalar 
Helmholtz equation in spherical coordinates. 

Another way to tackle this problem without consid- 
ering the vector harmonics employs the Hertz potential 
[il [l8j]. However, we prefer to adopt the same frame- 
work of Bohren and Huffman [6], in which the plane 
waves are directly expanded in terms of spherical vector 
harmonics. The solution of the scalar Helmholtz equa- 
tion is Vw {kr, cos 9, </>) = z n {kr)P™ (cos 8) exp(«m^), 
where z n (kr) is a generic Bessel spherical function and 
P™(cos#) are associated Legendre functions, n natural 
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and m integer. By means of ip nm , we can readily de- 
rive the spherical vector harmonics above defined. From 
the expansion of the incident fields in terms of M„ m and 
N nm , we find that only m — 1 contributes to this new 
representation due to the spherical symmetry of the scat- 
terer 0, 0] ■ Using the boundary conditions of this prob- 
lem, we can express the internal and the scattered fields 
in terms of spherical vector harmonics. In special, the 
coefficients of these expansions are referred to as the Mie 
coefficients. They provide the information about the in- 
teraction between the incident wave and the spherical 
particle. Explicitly, the boundary condition (r = a) is ex- 
pressed by (Ei+E s ~Ei) xe r = (Hi + H s -Hi)xe r = 0, 
where 1 is the index related to the particle (internal 
fields), and i and s refer to the incident and scattered 
fields, respectively; e r is the radial unity vector in the 
polar spherical coordinate system. 



A. Internal fields 

Assuming the incident EM wave is polarized in the 
direction, and the scattering center is placed at the origin 
of the coordinate system, we obtain an expression for the 
expansion of this field in terms of spherical harmonics. 
Imposing on the boundary between the sphere and the 
surrounding medium the continuity of the EM fields — 
in fact, their (electrical and magnetic) tangential compo- 
nents — , expressions are determined for the internal and 
scattered fields [a, 0, @1 . 

Using the same notation of Bohren and Huffman Q, 
we can give the components of the electric and magnetic 
vectors, Ei and Hi, respectively, of the interior field in 
a spherical coordinate system (r, 9, (f>) by 
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E n [d n Tn*pn{pi) ~ 7T„?/4 (/°1 )] 



factor) here adopted, the internal EM field (Ei,Hi) 
not the same presented in 



B. Magnetic Mie coefficients 

If we do not assume the equality between p and pi 
on the problem boundary condition, we can determine 
the magnetic Mie coefficients for the scattering (a„ and 
b n , obtained by Kerker at al.) and internal (c„ and d n ) 
fields [lHSj- Explicitly, 



d„ = 



fhlp n (mx)tjj n (x) - 1pn(xWn( mx ) 

fhij} n (mx)^ n (x) - £, n (xWn( mx ) 
ipn(mx)ip n (x) - fhip n (x)^' n (mx) 

i/j n (mx)^' n (x) - 7fi4 n (x)ip' n (rnx) 
mi 

ip n (mx)^' n (x) - m4 n (x)ip' n (mx) 
mi 

mip n (mx)£' n (x) - £, n {x)^' n {mx) 



(2) 

(3) 
(4) 
(5) 



with the assumption that the function domains are re- 
stricted in such a manner that the denominators do not 
vanish. The quantity x — ka is the size parameter of 
the spherical particle, being a its radius and k = |k| 
the wavenumber of incident and scattered waves, and 
£n(x) = x[j n {x) + iy n (x)] is the Ricatti-Hankel function 
of first kind. In addition, m = (p\t\/ pe) 1 ^ 2 is the rela- 
tive refraction index and fh ~ (pei/ P\t) is the relative 
impedance between the media. For p — pi, then rh = m 
and the usual expressions for the Mie coefficients ([5][5]) 
are recovered @, 0, [1] • 

There are some notation differences between this work 
and the one presented by Bott and Zdunkowski (l7j . 
Here, we use the same framework of Bohren and Huff- 
man Q, which have treated the scattering problem of 
light by means of International System of Units (SI), and 
have adopted exp(— iuf) as time-harmonic dependency 
for the fields. Otherwise, [I?} have used the same nota- 
tion as van de Hulst Q, which has dealt with the scat- 
tering problem in the Gaussian System of Units, and has 
adopted exp(icji) as time-harmonic dependency. These 
approach differences appear explicitly in the choice of 
the Hankel functions, which is strictly associated with 
the asymptotic limit for the scattered fields (the well- 
known far- field approximation), and consequently deter- 
mines the dependencies of the Mie coefficients. Another 
difference between these representations is related to the 
signal of the imaginary part of the relative refraction in- 
dex m = m r + imi, which is positive in the framework 
we have chosen [f| 0] ■ 



with p\ = k\r, E n — i n Eo{2n + l)/[n(n + 1)], n n and 
T n are angular functions defined in the appendix I A 21 and 
ipnipi) = Pijn(pi) is a Ricatti-Bessel function. The func- 
tions c n and d n are the internal magnetic Mie coefficients, 
which are presented in the section below. We outline 
that, because of the notation (units in the SI and time 



III. TIME-AVERAGED ELECTROMAGNETIC 
ENERGY 

For a linear, homogeneous, and isotropic medium, the 
classical theory of electromagnetism provides an expres- 
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sion for the time-averaged EM energy as an integral of 
the component intensities within the volume under anal- 
ysis. In the case of a spherical particle with radius a and 
internal complex EM field (Ei,Hi), we have the follow- 
ing relation |17l.ll9l]: 



W(a) 



xRe 



1 />a 

d(cos 9) I drr 2 
-l Jo 

(\E lr \ 2 + \E W \ 2 + \E M 



Mi 



(\Hi 



(6) 



In Eq. ©, with respect to the field representations, we 
outline that the permutation between a definite integral 
and a sum of an infinite series is not trivial. In the fol- 
lowing calculations, we are not concerned about showing 
explicitly each one of the simplifications; we just assume 
that the function series related to the field intensities 
converge uniformly in the domain 0<r<a,0<8<ir, 
< (f> < 2ir. Obviously, this mathematical condition is in 
agreement that the energy within a finite sphere is also 
finite. 



A. Electric and magnetic internal fields 

Looking closely to the definition of Eq. ([6]), one can 
ask about the contribution to the total internal energy 
associated with electric and magnetic fields separately, 
or even about the contribution of their fields compo- 
nents in spherical coordinates (r, 9, <j)) to this average en- 
ergy. These questions have not addressed by Bott and 
Zdunkowski in their (l7j paper. 



1. Radial component 

From Eq. ©, we obtain that the contribution of the 
radial component associated with the electric field to the 
internal energy is given by 

W rE (a) = / d(costf) &4> f \\rr 2 \E lr \ 2 

4 J-x Jo Jo 



a oo Co 



d„.dt 



Proceeding in the same way, one derives an analogous 
expression for the radial component related to the mag- 
netic field: 

W rH (a) = *-\Eo\ 2 ^t±Y.n{n + l) 

n=l 

2 / .1- f„ \|2 



xnn'(ra+l)(n' + l) / d(cos#) sin 2 6ir n n r , 



Eq. 



x K| 2 / dr\ Jn ( Pl )\< 
Jo 



x(2n+l)|c„|" / dr\ Jn ( Pl )\' . (8) 
Jo 

An important point to be noted here is that the integral 
above cannot be simplified by means of Eq. (|A1|) and 
recurrence relations presented in Appendix IA f [ 



2. Angular components 

Because of spherical symmetry of the system, it is not 
possible to write the contributions of the angular and 
azimuthal components to internal energy separately. If 
one tries to do that, the necessary relations to simply the 
double sums, as exemplified in Eq. (|7J|, do not appear. 
Fortunately, if one considers both (8, <f>) contributions to 
internal energy, these relations are not lost. Thus, using 
the relations (|A3|) and (|A4[) from Appendix IA 21 and the 
first term of Eq. ©, it follows that the time-averaged 
energy associated with angular components of the electric 
field is expressed by 

[W eE + W^ E ](a) = 5^ J' d{cos6) J\\$ 



x f drr 2 (\E ie \ 2 + \E^\ 2 
Jo 

= ?|£o| 2 ^X> + l) 



71=1 



x / dr(\c n Mpi)\ 2 + \dn^ n (pit 
Jo 



(9) 



Similarly, for the magnetic field we obtain 



[W eH + W^ H }(a) = -\Eq 



X)(2n+1) 



fci| 2 



(7) 



x / dr(KV™(pi)| 2 + |c„<(pi)| 2 ). 
Jo 

(10) 

Here, the same problem of the expression W r (a) arises. 
Whereas the integral of Ricatti-Bessel function is only 
another way to write Eq. (|A1[) . the second integral above 
cannot be simplified. 



B. Time-averaged internal energy 

From the expressions obtained in the previous section 
for each one of the internal field components, we can cal- 
culate the total time-averaged energy inside the sphere. 
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For the internal electric field, the expression is 



C. Dielectric sphere 



We (a) = W rE {a) + [W eE + W+ E ](a) 

= -W Re{mm)J2{(2n+l)\c n \ 2 l n (y) (11) 



+ \d n \ 2 [nl n+1 {y) + {n + l)X„-i(y)] } 



where 



1 f a 
Zn{y) = -o / drr 2 \j n 
a d Jo 



(Pi)| 



(12) 



is given by Eq. (IA1|) and Wo denotes the time-averaged 
EM energy of a sphere with radius a having the same EM 
properties of the surrounding medium: 



W 



'-na 3 \E \ 2 e 



(13) 



For sake of simplicity, the dependence of X n with respect 
to y* = m*ka, like the case of function W = W(a, y, y*), 
is omitted. 

In the same way, for the internal magnetic field, the 
average internal energy is given by 

W H (a) = W rH {a) + [W 8H + W^ H ]{a) 

= ^W Re(mm*) £ {(2n + \)\d n \ 2 X n {y) (14) 

n=l 

+ \c n \ 2 [nl n+1 (y) + (n+ l)l n -i(y)]} . 

Once we have expressions for electric and magnetic en- 
ergy within a sphere, it is possible to determine the ex- 
pression for the total time-averaged EM energy inside the 
scatterer: W(a) = We{o) + Wn(a). Explicitly, 



W(a) = -W J2\My)r 2 [nf3 n Z n+1 (y) 



where 



4 " (15) 

+ (n+ l) y n 2 n _i(y) + (2n + l)a n l n (y)] , 



\Mv)\ 2 [Re(mm)\c n \ 2 + Re (mm*) \d n \ 2 ] (16) 
\Mv)\ 2 [Re(mm)|d„| 2 +Re(mTO*)|c„| 2 ] (17) 



Also, to obtain analogous expressions to the ones pre- 
sented in |17| . Eq. ([TB]) can be rewritten as 



n=l y y 

A n (y*) A n {y) 



My*) My) 
y y* 



y y 

with y = mka and A n (y) = d y ]nip(y) 



(18) 



A particular situation to be considered here refers 
to a dielectric sphere, which has studied by Bott and 
Zdunkowski [l7|. With this aim, consider the non- 
magnetic case, i.e., (j, = /j,±. Thereby, it results that 
m = fh = (ei/e) 1 / 2 . With this assumption, note that 



Re (m 2 ) = (m 2 



2 )/2 and Re(r 



Substi 



tuting these into Eqs. and ((TTJ), the expression for 
the internal energy obtained by Bott and Zdunkowski is 
recovered. Once again, we emphasize that our notation 
is not the same that is employed in [l7j . Indeed, we can 
recover the same results by means of the following substi- 
tutions: ^ n (x) ( n (x), c n ---> md n , d n — » mc„, and 
assuming m = fh. Here, ( n (x) — x[j n (x) — iy n (x)] is the 
Ricatti-Hankel function of second kind, which is related 
to the choice of the time-harmonic dependence for the 
EM fields, like it is mentioned in the beginning of this 
description [(l-fioj. 

Employing the recurrence relations involving Bessel 
spherical functions 0, [2l[ , we obtain the derivative of 
first order A' n (y) = -I - A 2 n {y) + n(n + l)/y 2 . Therefore, 
using the L'Hospital rule, the limiting case of a perfect di- 
electric sphere, which takes place when m, — > 0, provides 
4y 2 lim mi ^ W(a) = 3W ET=i 7n(2n + 1)[1 + A 2 n (y) - 
n(n + l)/y 2 }, where 7„ = m 2 \^ n {y)\ 2 (\c n \ 2 + \d n \ 2 ). Un- 
less some commented notation differences, this result is 
the same obtained in (l7| . 



D. Absorption cross section 

The classical Mie theory provides a set of useful ex- 
pressions to calculate the scattering, total and absorp- 
tion cross-sections in the scattering process. Explicitly, 
call CT sca the scattering cross section and <7tot the extinc- 
tion (or total) cross section. Using the same framework 
of @, one can write 

2tt 



<r sca - ]T(2n + 1)(K| 2 + |&„| 2 ) , (19) 

(20) 



k 2 
2tt 



71=1 
OO 



otot = -p-^(2n+l)Re{a„ + 6„} . 

n=l 

Consequently, the absorption cross section cr a b s asso- 
ciated with the scatterer is defined in terms of cr sca and 
CTtot by the relation: cr abs = cr tot ~ cr sca - In other words, 
the absorption cross section in the Mie single scattering is 
determined by quantities and coefficients related only to 
the scattered EM fields @, 0, Although it is suitable 
and even natural to express the absorption cross section 
in terms of the internal coefficients c n and d n , notice that 
we do not do any reference to the internal EM fields. 

From the boundary conditions in the sphere problem 
Q , the Mie coefficients are linked by the equations below: 

h^\x)b n = j n (x) -j n (mx)cn , (21) 
h^\x)a n = j n (x) - fhj n {mx)d n . (22) 
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Thus, substituting the coefficients a n and b n into cr a b s = 
a tot— G sea an d manipulating that, we obtain the following 
expression respective to absorption cross section: 



Cabs 



n = l 



|V> ra (mx)|' 
|m£ n (a:)| 2 



ip n (mx) 
m&(x) 



\rndn 



(c„ + md n ) 



(23) 



Finally, using the definition of c n and d n given by 
Eqs. Q and © and the fact that Re[-i£*(x)£' n (x)\ = 
XnWV'nC 1 ) - ^(aOx^) = 1, where x„(x) = -xy n (x) 
is the Ricatti-Neumann function, we obtain 



2tt 



xlm 



(24) 



The exact expression (|24p for the magnetic absorption 
cross section in terms of the internal Mie coefficients is 
not found in the classical books of the scattering the- 
ory 0, 0, and it had not been determined so far. 



IV. NUMERICAL RESULTS 

In this section we present some numerical analysis from 
the exact equations determined in the sections above. 
Our aim is not to restrict our studies in some particu- 
lar case of magnetic scattering, but to introduce a gen- 
eral formulation of the internal energy which can be used 
whether in magnetic case {p ^ or in nonmagnetic 
one (jjl = iii). Here, all numerical results are obtained by 
means of a program created by us using the free software 
for scientific computation Scilab 5.1.1. As usual in the 
numerical Mie scattering, rather than to do infinite sums 
in Y]""^ in the calculation of the scattering quantities, 
which is impossible, we assume an approximation: finite 
sums with upper limit n max = x + 4a: 1 / 3 + 2, where x is 
the size parameter [8]. 



A. Normalized quantities 

For numerical studies, it is suitable to define dimen- 
sionless quantities related to the internal energy and the 
Mie coefficients: 

Wr m (m, fh, ka) = WE(m,m,kor,e,a) 

Wo(e, a) 

WT m (m, m, ka) = g^ffijgjlg) , (26) 

Wo{e, a) 

where We, Wh and W are expressed by Eqs. (fTTjl . ([T3|) 
and (fT3"|) , respectively. The dependence of We and Wh 
on the quantity m* is omitted. 



Therefore, one can define the normalization of the 
total time-averaged internal energy by the relation 

^norm = ^norm + J^norm Qr directly from Eq. (fT8]l : 

W^ m (m,fh,ka) = W(a)/W . 

Also, from [21j, we can obtain the recurrence rela- 
tion(2n+l) 2 | Jn ( Pl )| 2 = |pi| 2 {|j„-i(pi)| 2 + b„+i(pi)| 2 + 
2Re[j„_i(pi)j„ + i(p 1 )]}; thus, Eq. ([7]) can be rewritten as 



W rE (a) 3 



Wn 



= -Re(mm)^ 



n(n + 1) 
2n + l ' 



x |l„_ 1 (?/) +l n+1 (y) 
2 



(27) 



+ 



• / drr 2 Re[j n - 1 (pi)j n+1 (p* 1 )] 
Jo 



Note that the integral that appears in the sum above is 
quite similar to that one expressed in Eq. (|A1|) . Although 
this one cannot be simplified by means of Eq. (|A1|) . it 
is possible to show numerically that the result of this 
integral is proportional to a 3 . It means that we can study 
the radial contribution to the internal energy normalized 
by Wo using only the dimensionless parameters m, fh and 
ka. The same argument can be applied to both [WgE + 
W < j > E)(a)/Wo and the analogous expressions respective to 
the internal magnetic field, given by Eqs. (JSJ) and (|T0|) . 

In the situations considered here, the values of Re(mm) 
and Re (mm*) are very close in such a way that We(ol) ~ 
Wh(cl). Thus, we only consider the total time-averaged 
internal energy in our analysis. Further, we remark that 
although in these studies the relative magnetic perme- 
ability is assumed to be real, there is no such restriction 
in the calculated expressions. For soft microwave ferrites, 
a more realistic study should consider the magnetic loss. 

Figs. [1] and [2] illustrate a succession of narrower picks 
of the values of the normalized internal energy in single 
magnetic Mie scattering as a function of the size param- 
eter and relative magnetic permeability. Here, we are 
not concerned about to study in details the resonances 
picks and ripple structure due to the internal coefficients 
c„ and d n \TA I22I [23j . The internal energy of a mag- 
netic sphere presents resonances peaks even in the limit 
of small geometric size (compared to the wavelength). 
These resonances are due to the increase of the total 
cross section due to magnetism. This increase leads to a 
decrease of the photon mean free path in multiple scat- 
tering regime in a disordered system. The smaller mean 
freepath favors the localization phenomenon as pointed 
in [3j-l5j. For dielectric spheres, these narrower resonance 
picks are well known and they are referred to as mor- 
phology dependent-resonances (MDRs) In the Mie 
theory, for large size parameters, these MDRs are com- 
monly observed at the scattered and internal intensities 
and at the total cross section. 

In our system, observe that when one increases the 
contribution of the magnetism in the scatterer, the values 
of the internal energy W{a) become much larger than 
Wq, and the narrower picks appear even for small size 
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Figure 1: Comparison among the distributions of W(a)/Wo 
in function of the size parameter ka. The values of the rela- 
tive permittivity and permeability related to a non- absorptive 
sphere are ei/e = 1.4161 and = 1, 10, 100, 1000, respec- 

tively. The internal energy W(a)/Wo is calculated in the in- 
terval 10" 6 < ka < 1, S(ka) = 10" 4 . 
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Figure 2: The normalized time-averaged internal energy 
W(a)/Wo within a non-absorptive magnetic sphere plotted 
as a function of the size parameter ka. The values of the 
relative permittivity and permeability are ei/e = 1.4161 and 
Hi/fi = 10 4 , respectively. The internal energy W(a)/Wo is 
calculated in the interval 10 -6 < ka < 1, S(ka) = 10 -4 . 



geometrically the scatterer is much smaller compared to 
the wavelength. In other words, the incident EM wave 
interacts strongly with the optical cross section instead 
of the geometrical one 



10 b 
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Figure 3: The normalized time-averaged internal energy 
W(a)/Wo plotted as a function of the size parameter ka. The 
values of the relative permeability are = 1 (nonmagnetic 
sphere) and /ii/p = 100 (magnetic sphere). The relative re- 
fraction index is m = 1.334+1.5 x 10 _9 i, which have been used 
in [ItJ in the nonmagnetic scattering approach. The quanti- 
ties are calculated in the interval 1 < ka < 50, 8(ka) = 0.01. 



It must be mentioned that this assumption of a non- 
absorptive magnetic diffusor is conditioned to the fre- 
quency range of the incident beam. Usually, it is con- 
trolled with an external static magnetic field |llT - il3l [l5| . 
Indeed, there is a wide variety of soft ferrites which ex- 
hibits very large values of relative magnetic permeability 
at applied fre quen cies typically below 100 MHz with low 
magnetic loss [TJ. For sake of simplicity and generality, 
the situations considered here do not take into account 
the scatterer in the magnetized state, and a scalar value 
for the [ii/n is adopted. The dependence on the angular 
frequency lj of the incident EM wave remains implicit on 
the value of the size parameter ka. Given a value of oj 
for the incident EM wave, a surrounding medium (e, /i) 
and a scatter (ei,/xx), one readily obtains k = uj(iie) 1 / 2 
and k\ = mk (see Fig. 0]). 



B. Weak absorption regime 



parameters. In the nonmagnetic case reported in [17j . 
there is an oppositive tendency, that is, both the internal 
energy and the absorption efficiency enhance with the 
size parameter. These difference between a nonmagnetic 
case and a magnetic one is illustrated in Fig. G3 This is 
due to the increase of the total cross section even though 



In the weak absorption regime (wa), it is quite evi- 
dent the relation of the time-averaged internal energy 
and the absorption efficiency when we compare them, as 
it is shown in the Fig. [5l This correlation between these 
quantities in nonmagnetic scattering have been studied 
in [13]. 

Analytically, for mi m r and rhi m r , we can write 
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Figure 4: The separation of the total time-averaged internal 
energy in radial and angular contributions respective to both 
electric and magnetic fields. The values of the relative elec- 
tric permittivity and magnetic permeability are ei/e = 10 
and jUi/p = 100, respectively. The component contribu- 
tions W r (a)/Wo and We,<p(a)/Wo are calculated in the in- 
terval 10~ 6 < ka < 2, 5(ka) = 10" 4 . 

yi _y*i ^ ^ lx 2 m r mi, Re(mm) sa m r fh r and Re (mm*) « 
m r m r . Using these approximations in Eq. (|18p. it follows 
that 

/ n 3 TT , m r 2fh r ,„ 
W wa (a) « -W _!:_JL^(2n-|-l 
8 m ? - rm: ^— ' 

' n— 1 

x (Ic^ + ldJ^Im^W^y*)] -(28) 

Once the absorption efficiency in the Mie single scattering 
is defined by Q a bs = Cabs/cg, where <r g = ira 2 is the 
geometric cross section and <7 a b s is expressed in Eq. (1241) , 
we can write that 

W wa (a) « |w — Qabs , (29) 
8 m^x 

which is the same relation obtained in [l7j in the nonmag- 
netic case. Indeed, this approximation is valid wherever 
m, <C m r and m 2 ; <C m r ; it is not affected by the value 

of fil/fl. 

In addition, for the case in which m r s=s 1, one obtains 
W / wa (a) « Wo- Therefore, in this particular situation, 
one can write Q a bs ~ 8x/(3m»), which is a well-known 
expression 0, fl7j . 




Figure 5: Ratio between the absorption efficiency Q a b s and 
the normalized time-averaged internal energy W(a)/Wo plot- 
ted as a function of the size parameter ka. The values of the 
relative permeability and refraction index are /Lti//i = 1 (non- 
magnetic sphere) and m = 1.334 + 1.5 x 10~ 9 i, respectively. 
The quantities are calculated in the interval 1 < ka < 49, 
5(ka) = 2. The angular coefficient of linear regression is ap- 
proximately 2.997 x 10 -9 , which is in agreement with Eq. (I29|) : 
8mi/(3m r ) « 2.998 x 10" 9 . 



expressions for the contributions to the EM internal en- 
ergy related to the fields components separately. The ex- 
pressions for the EM internal energy within a dielectric 
sphere and the relation derived in |17| in the weak absorp- 
tion regime between the internal energy and the absorp- 
tion efficiency are recovered. In special, we find that the 
magnetism of the particle does not break the linear re- 
lation between the absorption efficiency and size param- 
eter. To do so, we analytically calculate an expression 
for the absorption efficiency, which depends only on the 
internal magnetic Mie coefficients. In addition, we cal- 
culate the limiting cases of the magnetic Mie coefficients 
and present some important properties of the radial func- 
tions which are used to simply the obtained expressions. 
Finally, the main result of this work is that, even for 
small scatterers compared to the wavelength, the value 
of the EM internal energy within a magnetic sphere is 
much larger than that one associated with a sphere with 
the same properties of the surrounding medium. Physi- 
cally, we ascribe this fact to the enhancement of the total 
cross section due to the magnetism in the scatterer. 



CONCLUSION 



In this paper we generalize the exact expression of 
the time-averaged EM internal energy, obtained firstly 
in [T3|, to the case of magnetic spherical scatterers. Us- 
ing the same framework of (| and assuming the mag- 
netic scattering approach [l|, we determine analytical 
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Appendix A: Special functions 
1. Radial functions 



we imperatively have to consider n — 1 and n 
keep consistent orders in the Mie coefficients. 



2 to 



For the situation in which rrii =^ is verified, that is, 
the imaginary part of relative refraction index (absorp- 
tive component) is not null, we can write 



?" 2 b'n(pi)| dr 



a 3 [y*3n(y)j' n (y*) ~ yjn(y*)j' n (y)] 



i*2 



2a 3 \j n (y)\ Re 



- y 

<Pn(V*) 



y - y 



(Al) 



where pi = mkr, y = mka, ip n (y) = yd y [ln^ n (y)]. 
Eq. (|Alj) is provided, in terms of Bessel cylindrical func- 
tions, by Watson [20| . In the present context, to treat 
only with spherical Bessel functions, we have used the 
relation (2p 1 /n) 1 / 2 j n (p 1 ) = J n+1 , 2 (pi) M- 

Also, if the relative refraction index m is real, accord- 
ingly to [2(J, the integral in Eq. (jAll) can be simply 
rewritten as 

a 3 

r 2 ]l(fii)Ar = — [j 2 n (y) - j n -i(y)jn+i(y)] , (A2) 

which is obtained by taking the limiting mj —> in 
Eq. (|A1[) . and using L 'Hospital rule and recurrence re- 
lations. 



2. Angular functions 

In the expansion of EM fields, it becomes natural to 
define the angular functions 7r„(cos#) = P* (cos 9) /sin 9 
and r n (cos#) = dP*(cos0)/d0, where 9 is the scattering 
angle and P^ is an associated Legendre function of first 
order. These angular functions are quite convenient in 
the calculation of fields intensities. 

Due to properties involving the associated Legendre 
functions, 7r„ and r„ satisfy the following expressions, 
Vra, n' E IN: 



d^osfl) (n n Tr nl + T n T n <) = 



2n 2 (n + l) 2 
2n + 1 



d(cos6<) (7r„r„/ 



d(cos ( 



2n(n + 1) 



sin 2 9 



2n+ 1 



(A4) 



.(A5) 



These expressions facilitate the determination of quan- 
tities involving fields intensities. We emphasize that, in 
the classical books of scattering _theory, Eq. (IA5I) is not 
found in this explicit form [T 



1. Small particle limit 



For the small argument limit into the Mie scattering 
coefficients, we obtain 



a i 



tx 3 ip i (mi) — 2mm 
3 (pi(mx) + mm 
ix b [(fi (mi) — mm 
5 



mmip\(mx) 



[ipi (mx) + mm] 



x 

~~9 



ifi(mx) — 2mm 



tpi(mx) + mm 
2m/rh 



0(x 7 ) 



(Bl) 



ix tpi(mx) 
3 Lpiijnx) + m/fh 



ix 



[(pi(mx) — m/m] 2 — (m/m)(p 1 (mx) 



fi r 



[ipi(mx) 
(pi(mx) — 2m/fh 



m/m] 



a 2 
b 2 



ipi(mx) + m/m 
ix 5 if2 (mx) — 3mm 
45 if2 (mx) + 2mm, 
ix 5 ip 2 (mx) — im/rh 
45 Lp2(mx) + 2m/fh 



+ 0(x 7 ) 
-0(x 7 ) , 
+ 0(x 7 ) . 



(B2) 
(B3) 
(B4) 



For the Mie internal coefficients, the approximations 
assume the form below: 



ci 



mx 



m/m 



(-2 



ipi(mx) [tpi(mx) + m/m] 

mx 4 (m/m) [ipi(mx) — m/m 
ipi(mx) 2[ipi (mx) +m/fh] 2 
(m/fh)x 2 mm, 
ipi(mx) [ipi (mx) + mm,] 
(m/m)x A mm, [ifi(mx) — mm] 
ipi(mx) 2 [ipi (mx) + mm] 2 
m (m/m)x 3 
3f/'2(wx) [if2(mx) + 2m/rh] 



(B5) 
+ 0(x 5 lV6) 
(B7) 
0(x 5 ) , (B8) 
+ 0(x 5 ) , (B9) 



(m/m) 



mmx 



3ip2 (mx) [ip2(mx) + 2mm\ 



+ 0(x 5 ) 



(BIO) 



Appendix B: Limiting cases 

In these particular cases, we remark that n = 1 is suf- 
ficient to study the nonmagnetic scattering theory. Here, 



Note that, for these approximations, the scattering coef- 
ficients ai, d2, bi and 62 have order O \(ka) 7 ] , whereas 
the internal coefficients c l5 C2, c?i and d 2 are O [(ka) 5 ]. 
Terms for n > 2 are ignored here. 
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2. Rayleigh approximation 

In this approximation, in which \m\x <C 1, the Mie 
scattering coefficients can be write as 



«i 



-2ix mm — 1 



3 mm + 2 

ix 5 m 3 m — 6mm, + (mm) 2 + 4 
5 (mm + 2) 2 

9 \ mm + 2 J 
-2ix 3 m/m — 1 



(BH) 



3 m/m + 2 

?a; 5 m 3 m — 6m/m + (m/m) 2 + 4 



4x 6 / m/m — 1 x 



«2 
&2 



9 \m/m + 2 
-zx 5 mm, — 1 



(m/m + 2) 2 

+ 0(x 7 ) , 



15 2mm + 3 

-IX 5 m/m — 1 



+ 0(x 7 ) , 



15 2m/m + 3 
and the Mie internal coefficients assume the form 



(B12) 
(B13) 
(B14) 



ci 
di 



2m + m 
_3x^ 

2 

3 



1 



(mx)' 



1 



10 

(mi) 2 



2 + mm 
_3x^ 
2 



1 



10 

(mx) 



(2m — m) 
(2m + m) 2 



1 + 



10 

(mx) 2 



10 



mm(3 + 2m/m) 



(2 — mm) 
(2 + mm) 2 

0(x 5 ) , 



<D(x 5 )(pi5) 



0(x 5 )(pi6) 



m(3 + 2mm) 



C(x 5 



(B17) 
(B18) 



Taking the particular case m = m, Mie coefficients for 
nonmagnetic scattering are recovered @. 



compared to wavelength. Using the expressions for large 
argument limit present in [25] ]. one can obtain 



a i 



tx 3 xtan(mx) + 2m 
3 xtan(ma;) — m 

ix 5 [xtan(mx) + m,] 2 + mxta,n(mx) 
5 



[x tan(mx) — m] x 



x 

~~9~ 



x tan(mx) + 2m 



xta,n(mx) — m 
ix 3 mx tan(mx) + 2 
3 mx tan(ma;) — 1 



+ 0(x 7 ) , (B19) 



ix 5 [mx t&n(mx) + l] 2 + mx ta,n(mx) 
~5 



x 

~~9 



mx tan(mx) 



mxtan(ma;) — 1] 

+ 0(x 7 ) 



mx tan(mir) — 1 



«2 



tx 5 x — 3mtan(mx) 
45 x + 2m tan(mx) 
ix 5 mx — 3 tan(mir) 
45 fhx + 2 tan(mx) 



the internal coefficients are 



ci 



+ 0(x 7 ) , 
+ 0(x 7 ) ; 



(B20) 
(B21) 
(B22) 



di 



C2 



cos(mx) [mx tan(mir) — 1] 
mi 4 [mxtan(ma;) + 1] 
2cos(mx) [mxtan(ma;) - l] 2 
m x 2 
cos(mx) [xtan(ma;) — m] 

mi 4 [x tan(mx) + m] 
2cos(mx) [xt&n(mx) -fhf 
—m x 3 
3 cos(mx) [2 tan(mx) + rni] 

— m x 3 
3cos(ma;) [2mtan(mx) + .t] 



0(a; 5 )(B23) 



0(x 5 ) ,(B24) 



C(x 5 ) , (B25) 



0(x 5 ) . (B26) 



3. Ferromagnetic limit for x -C 1 

This approximation, similar to Rayleigh limit, is de- 
rived directly from the approximation of small spheres 



In this case, since low order in the size parameter is 
used, one can obtain an analytical expression for the 
physical quantities, such as cross sections, for instance. 
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